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Introduction: The quest for 
protected qubits 



The basic building block of quantum computation is the qubit, a system with two 
(nearly) degenerate states that can be used to encode quantum information. Real 
systems typically have a full spectrum of excitations that are considered illegal from 
the point of view of a computation, and lead to decoherence if they couple too strongly 
into the qubit states during some process. See Fig. 11.11 The essential problem then 
is to preserve the quantum state of the qubit as long as possible to allow time for 
computations to take place. 

Assuming the gap A to the illegal states is reasonable, we can quite generally 
describe the dynamics of the qubit state by an effective Schodinger equation 



where H^g is the effective qubit Hamiltonian. In quantum optics. Hcs is often known 
with high precision. This is not so in condensed matter systems such as quantum dots. 
Even worse, fleff may fluctuate or include interaction with the environment. This 
causes decoherence of the qubit state. 

Ideally, we would like to arrange for H^s to be zero (or Hes = £l) for some good 
reason. Usually, we use a symmetry to protect degeneracies in quantum systems. For 
example, a quantum spin i has a two-fold degeneracy protected by the SU (2) symme- 
try, as do the 2s + 1 degeneracies of higher spins s. Indeed, any non-Abelian symmetry 



^|vE') = -*i7eff|*) 



(1.1) 
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Fig. 1.1 Spectrum of a physical qubit system. 
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would work. Unfortunately, the SU{2) symmetry of a spin is lifted by magnetic fields 
and it's generally difficult to get rid of stray fields. 

Rather than symmetry, in what follows we will look to topology to provide us 
with physically protected degeneracies in quantum systems. In particular, we will 
examine a number of exactly solvable models in one and two dimensions which exhibit 
topologically phases - that is, gapped phases with a protected ground state degeneracy 
dependent on the topology of the manifold in which the quantum model is embedded. 
In Sec. [2] we warm up with the study of several quantum chains that exhibit Majorana 
edge modes and thus a two-fold degeneracy on open chains. The topological phenomena 
available in two dimensional models are much richer and will be the focus of the 
remaining three sections. We introduce and solve the toric code on the square lattice in 
Sec. [31 exhibiting its topolological degeneracy and excitation spectrum explicitly. The 
following section steps back to examine the general phenomenology of quasiparticle 
statistics braiding in two dimensional models. Finally, in Sec. [5] wc introduce the 
honeycomb lattice model which exhibits several kinds of topological phases, including 
that of the simple toric code and, in the presence of time reversal symmetry breaking, 
a gapped phase with chiral edge modes protected by the topology of the Fermi surface. 
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Topological phenomena in ID: 
boundary modes in the Majorana 
chain 



We will consider two examples of ID models with Z2 symmetry and topological 
degeneracy: the transverse field Ising model (TFIM) and the spin-polarized supercon- 
ductor (SPSC). Although these models look rather different physically, we will find 
that they are mathematically equivalent and that they both exhibit a topological phase 
in which the ground state degeneracy is dependent on the boundary conditions of the 
chain. That is, the ground state on an open chain is 2-fold degenerate due to the pres- 
ence of boundary zero modes, whereas the ground state is unique on a closed loop. 
This topological degeneracy will be stable to small local perturbations that respect 



the Z2 symmetry. More details on these models may be found in Kitaev (2000 1 . 

1. The transverse field Ising model is a spin-1/2 model with Hamiltonian: 

N-l N 

Hs = -jY.^l^U-h^T.-r (2-1) 

Here J is the ferromagnetic exchange coupling in the x direction and is a 
uniform transverse (z) field. This model has a Z2 symmetry given by a global 
spin flip in the a-x basis: 

N 

Ps = X{-' (2.2) 

2. The spin-polarized 1-D superconductor is a fermionic system with Hamiltonian: 

N-l 

Hp = (^-?i;(a]aj+i -|- a]_|_j^aj) + Aaja-j+i + A*aj^^af^ 

j=i \ / 

where aj and are fermionic annihilation and creation operators, w is the hop- 
ping amplitude, A is the superconducting gap and /z is the chemical potential. 
For simplicity, we will assume that A = A* = so that 
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N-l N 



Hf^wY, {a, - a]){a,+i + a]+i) - - I/2) • (2.4) 

This model has a Z2 symmetry given by the fermionic parity operator: 

Pf = (-1)^^"^*"^' (2.5) 

Although the two models are mathematically equivalent, as we will sec in Sec. 
they are clearly physically different. In particular, for the superconductor, the Z2 
symmetry of fermionic parity cannot be lifted by any local physical operator, as such 
operators must contain an even number of fermion operators. Unfortunately, for the 
spin system the degeneracy is lifted by a simple longitudinal magnetic field J^j '^j 
and thus the topological phase of the TFIM would be much harder to find in nature. 

2.1 Nature of topological degeneracy (spin language) 



Consider the transverse field Ising model of Eq. (|2.ip . With no applied field, there are 
a pair of Ising ground states [h^ = 0): 

|V^^) = | — •••-), |^^) = | — •••-). (2.6) 

The introduction of a small field hz allows the spins to flip in the basis. In particular, 
tunneling between the two classical ground states arises via a soliton (domain-wall) 
propagating from one side of the system to the other: 

H^^...^)^|^:^^...^)^|^^:^...^) (2.7) 
^ I ^^^: . . . ^ ) ^ . I . . . ^ ). (2.8) 

As usual, the tunneling amplitude t associated with this transition falls off exponen- 
tially in the distance the soliton must propagate 

t ^ e-^/« (2.9) 

where ^ is the correlation length of the model. The two-fold degeneracy is therefore 
lifted by the effective Hamiltonian: 

H.. = (2.10) 

The splitting is exponentially small in the system size and the two-fold degeneracy is 
recovered in the thermodynamic limit as expected. Moreover, it is clear why introduc- 
tion of a longitudinal field will fully split the degeneracy. 

2.2 Reduction of TFIM to SPSC by the Jordan-Wigner 
transformation 



To show the equivalence of the one dimensional models introduced above, we will use 
a standard Jordan-Wigner transformation to convert the spins of the Ising model into 
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fermions. It is perhaps not surprising that a fermionic description exists for spin 1/2 
systems - we simplify identify the up and down state of each spin with the presence 
or absence of a fermion. The only difficulty arises in arranging the transformation 
so that the appropriate (anti)-commutation relations hold in each description. The 
Jordan- Wigner transformation docs this by introducing string-like fermion operators 
that work out quite nicely in 1-D nearest neighbor models. 
To reduce Hs to Hp, we 

1. Associate the projection onto the z-axis of the spin with the fermionic occupation 
number: 

|T)^« = 0, U)^n=l. (2.11) 

That is, 

cr| = (-1)°^^-. (2.12) 

2. Introduce the string- like annihilation and creation operators 




^7 (2.13) 

\k=l ) 

where and a~ are the usual spin raising and lower operators. At this stage, we 



can check that the usual fermionic anticommutation relations hold for the , a| : 

{a„a]}=% (2.14) 

3. Observe that 

so Hs (Eq. (EHI)) reduces to Hp (Eq. with 

w = J, ^ = -2/i^ (2.16) 

2.3 Majorana operators 

Majorana operators provide a convenient alternative representation of Fermi systems 
when the number of particles is only conserved modulo 2, as in a superconductor. 
Given a set of TV Dirac fermions with annihilation/creation operators aj,aj, we can 
define a set of 2TV real Majorana fermion operators as follows: 

C2j-i = flj + a] 
a," — aJ 

C2, = (2.17) 



These operators are Hermitian and satisfy a fermionic anticommutation relation: 
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Fig. 2.1 Majorana chain representation of 1-d superconductor. Each boxed pair of Majoranas 
corresponds to one site of the original fermionic chain. 



Or, more compactly, 



4 = 1, CfcQ =~QCfe(fc^O- (2.18) 



{ck,ci} = 2Ski. (2.19) 



From any pair of Majorana operators, we can construct an annihilation and creation 
operator for a standard Dirac fermion (a = (ci +ic2)/2 and h.c), and thus the unique 
irreducible representation for the pair is a 2-dimensional Hilbert space which is either 
occupied or unoccupied by the a fermion. 

Both models Hs and Hp can be written as 

. / N N-1 
i^maj = ^ XI '=2j-lC2i + W ^ C2jC2j + l I (2.20) 



J = l i=i 



where v = ^ ^/^/2 and w = J. The Z2 symmetry of fermionic parity is given in 
the Majorana language by 



N 

^maj = n(-*^2fc-lC2fe). (2.21) 

fc=l 

We can view this model graphically as a chain of coupled Majorona modes, two to 
each of the N sites of the original problem as in Fig. 12.11 If v = 0, then the Majorana 
modes at the ends of the chain are not coupled to anything. This immediately allows 
us to identify the 2-fold ground state degeneracy in i/maj as the tensor factor given 
by the 2-dimensional representation of the boundary pair ci, C2n- 

We will see in Sec. 12.41 that ii v but |?;| < w, the operators ci and C2n are 
replaced by some boundary mode operators bi, br- The effective Hamiltoirian for this 
piece of the system is then 

-ffcff = ^ebibr = eia^a - i) (2.22) 

where e ~ e~^/^ and a, are the Dirac fermion operators constructed from the 
boundary pair. Thus, the ground state degeneracy is lifted by only an exponentially 
small splitting in system size. 
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2.4 General properties of quadratic fermionic Hamiltonians 

We now step back and consider a generic quadratic fermionic Hamiltonian: 



H{A) 



4 ^ 



-^j k Cfc 



(2.23) 



where ^ is a real, skew-symmetric matrix and the Cj are Majorana fcrmion operators. 
The normalization 



^ is convenient because it has the property that 
[-iH{A),-iHiB)] = ^tH {[A, B]) 



where A,B G 5o{2N), and H{A),H{B) act on the Fock space = 
provides a natural representation of so(27V). 
We now bring H{A) to a canonical form: 



(2.24) 
Thus H{-) 



i^canonical = ^ X! ^kb'kb'^ = ^ efc(al«fe ^ 2^ 



(2.25) 



fc=l 



k=l 



where b'j., b'^ are appropriate real linear combinations of the original Cj satisfying the 
same Majorana fermion commutation relations and the a^, a^. are the annihilation and 
creation operators associated to the foj., b'^ pair of Majoranas. This form for H follows 
immediately from the standard block diagonalization of real skew symmetric matrices 



A = Q 



/ ei 
-ei 



V 



£2 
-£2 
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Q^, Q eOi2N),ek>0 



(2.26) 



From this form it is easy to check that the eigenvalues of A are ±iek and that the 
eigenvectors are the coefficients of c in 5^, a^. 

If some of the vanish, then we refer to the associated fermions as zero modes. In 
particular, these will lead to ground state degeneracies, since occupation or nonoccu- 
pation of such modes does not affect the energy. For the Majorana chain of Eq. (|2.20[) . 
we have 

/ V \ 

-V w 
~w V 
—V w 
~w 

V •7 

We can find a vector u such that uA = by inspection 

2 



(2.27) 



V / V \ ^ 

w = (l,0,-,0, - ,o,---) 



(2.28) 



This vector leads to a left boundary mode 
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bi=^UkCk (2.29) 

while an analogous calculation starting at the right end will find a right boundary mode 
br- These modes form a Majorana canonical pair, leading to a two- fold degeneracy of 
the ground state of the chain. Clearly, Uk e~^/^ falls off exponentially from the edges 
of the chain with correlation length = In 1^1, as expected in section [2?3l 



2.5 Why are the boundary modes robust? 

In the simple case of a quadratic fermion Hamiltonian, we know that the modes cor- 
respond to eigenvalues of a skew-symmetric real matrix. These come in pairs ±ie, in 
general, and the case e = is special. In particular, if the pair of Majoranas corre- 
sponding to a zero mode are physically well separated, we expect perturbations to 
have trouble lifting the boundary degeneracy. 

More generally, for interacting fermions, we can extend the symmetry group Z2, 
generated by P = Pmaj, to a non-comm,uting algebra acting on the the ground state 
space £. First, in the noninteracting limit, at w = 0, we define 

3 

X = Y\[{-ic2k-iC2k) (2.30) 

fc=i 

where Y = C2j+i is a local Majorana operator at site 2j -I- 1. A straightforward calcu- 
lation shows that 

XP=-PX (2.31) 

and that [-ff, X] = so that the algebra generated by X,P acts on £ nontrivially. 
We now allow Y to vary as we adiabatically turn on interactions and, so long as an 
energy gap is maintained, we expect Y to remain a local operator near 2j, which we 
can separate from the boundary by suitably large choice of j. That is, to find 1", one 
needs to know the ground state or at least the structure of the ground state near 2j. 



This is a nontrivial operation but see Hastings and Wen (2005 ) for more details 



3 

The two-dimensional toric code 



The toric code is an exactly solvable spin 1/2 model on the square lattice. It exhibits a 
ground state degeneracy of 4^ when embedded on a surface of genus g and a quasiparti- 
cle spectrum with both bosonic and fermionic sectors. Although we will not introduce 
it as such, the model can be viewed as an Ising gauge theory at a particularly simple 
point in parameter space (sec Sec. 14. 5p . Many of the topological features of the toric 
code model were essentially understood by Read and Chakraborty (1989), but they 



did not propose an exactly solved model. A more detailed exposition of the toric code 



may be found in Kitaev (2003 ) 



We consider a square lattice, possibly embedded into a nontrivial surface such as 
a torus, and place spins on the edges, as in Fig. 13.11 The Hamiltonian is given by 



Ht — —Je ^ As — Jrn ^ Bp 



(3.1) 



where s runs over the vertices (stars) of the lattice and p runs over the plaquettes. 
The star operator acts on the four spins surrounding a vertex s. 
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Fig. 3.1 A piece of the toric code. The spins live on the edges of the square lattice. The 
spins adjacent to a star operator As and a plaquette operator Bp are shown. 
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As= II <y] (3.2) 

while the plaquette operator acts on the four spins surrounding a plaquette, 

Bp = n "j- (3.3) 

jedp 

Clearly, the Ag all commute with one another, as do the Bp. Slightly less trivially, 

A,Bp = BpA, (3.4) 

because any given star and plaquette share an even number of edges (either none or 
two) and therefore the minus signs arising from the commutation of and on 
those edges cancel. Since all of the terms of Ht commute, we expect to be able to 
solve it term by term. 

In particular, we will solve Ht working in the basis. Define classical variables 
Sj = ±1 to label the basis states. For each classical spin configuration {s}, we can 
define the plaquette flux 

"'p(s) = Yl (3-5) 
jedp 

If u>p — —1, we say that there is a vortex on plaquette p. 

3.1 Ground states 

To find the ground states |^) of Ht, we need to minimize the energy, which means 
maximize the energy of each of the A^ and Bp terms. The plaquette terms provide the 
condition 

Bpl^f) = 1*) (3.6) 

which holds if and only if 

l*)= E ^^l^> (3.7) 

{s:Wp{s) = l Vp} 

. That is, the ground state contains no vortices. The group of star operators act on 
the configurations s by flipping spins. Thus, the star conditions 

A,\^) =. I*) (3.8) 

hold if and only if all of the Cg are equal for each orbit of the action of star operators. 
In particular, if the spin flips of As are ergodic, as they are on the plane, all Cs must 
be equal and the ground state is uniquely determined. 

On the torus, the star operators preserve the cohomology class of a vortex-free spin 
configuration. In more physical terms, we can define conserved numbers given by the 
Wilson loop like functions 

Wi(s)=]^Sj, l^h^h (3.9) 

where h and I2 are two independent non-trivial cycles on the square lattice wrapping 
the torus (Fig. 13. 2p . Any given star will overlap with a loop / in either zero or two 
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Fig. 3.2 Large cycles on the torus. 
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Fig. 3.3 Electric and magnetic path operators. 

edges and therefore As preserves wi. Since there are two independent loops on the 
torus, each of which can have wi = ±1, there is a four-fold degenerate ground state: 



i*)= E 



\s). 



(3.10) 



{s:ujp(s) = l Vp} 



3.2 Excitations 

The excitations of the toric code come in two varieties: the electric charges and mag- 
netic vortices of a Z2 gauge theory. We will see this connection more explicitly later. 
In the following, we restrict attention to the planar system for simplicity. 
To find the electric charges, let us define the electric path operator 

w^/^)=n^i (3.11) 

where I is a path in the lattice going from si to S2 (see Fig. l3.3p . This operator clearly 
commutes with the plaquette operators Bp and with all of the star operators As except 
for at the end points si and S2, where only one edge overlaps between the star and 
the path and we have 

W^^'^As.^-AsM"^. (3.12) 

Therefore, the state 

l*.i....)-W-/'^l*o), (3.13) 

where |^'o) is the planar ground state, is an eigenstate of the Hamiltonian with excita- 
tions (charges) at si and S2 that each cost energy 2 Jg to create relative to the ground 
state. 
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An analogous construction will find the magnetic vortices: we can define a dual 
path operator 

w^/r' = n ^1 (3-14) 

where the path I* lies in the dual lattice (see Fig. 13. 3p and goes from pi to p2- In this 
case, the stars Ag all commute with as do all of the plaquette operators Bp 

except the two at the end points of I*, which anticomniutc. Thus, the W^l"'^ operator 
creates a pair of magnetic vortices on the plaqucttes pi and p2 at an energy of 2 J™ 
each. 



4 

Abelian anyons and quasiparticle 
statistics 



Let us discuss what can possibly happen if we exchange two particles in two dimensions. 
To ensure that particle statistics is well-defined, we assume that there is no long-range 
interaction and that the phase is gapped. If we drag two particles around one another 
adiabatically, 



Top View 




t 



Side View 



then we expect both dynamical phase accumulation and a statistical effect due to the 
exchange. We are well acquainted with this effect for everyday bosons and fermions, 
for which: 

Bosons: I'S) ^ 1*) 
Fermions: I'i') — I'J') 



(4.1) 



where we have dropped the dynamical phase so as to focus on the statistics. In both 
of these standard cases, a full rotation (two exchanges), 




Top View 




leaves Ivp) unchanged. 
In principle, 



ab 




(4.2) 



is an arbitrary phase factor or even an operator [braiding operator). If the two particles 
are distinguishable (a 7^ 6), then Rah does not have an invariant meaning, but the 
mutual statistics 
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Rba ■ R. 




ab 



(4.3) 
does. 

Let us illustrate this in the toric code model. In section [321 we found two kinds of 
quasiparticle excitations in the toric code: electric charges (e) and magnetic vortices 
(m). Since path operators of the same type commute with one another, it is easy to 
show that each of these are bosons. However, they have nontrivial mutual statistics. 

To calculate the mutual statistics, consider taking a charge e around a vortex m. 
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Let 10 be some state containing a magnetic vortex at pi. Under the full braiding 
operation, 

10- m-lj le) 



n 

inside / 



(4.4) 



where the second line is a Stokes' theorem like result relating the product around a 
loop to the products of internal loops. Since 



for the plaquette pi containing the vortex, we have that 

10- -10, 

or 

e m 



m 



(4.5) 
(4.6) 




(4.7) 
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Using the bosonic self-statistics equations, 




e 



e 




m m 



ee ee m m m m (4.8) 
we can derive the nontrivial corollary that composite e — m particles are fermions: 



e m 



e m 




em em 



e m 



e m 




m 



e m 



m 



e m 



m 



(4.9) 



4.1 Superselection sectors and fusion rules 



Initially, we exhibited two kinds of bosonic excitations in the toric code model (charges 
e and vortices m) in the solution of the Hamiltonian. After a bit of work, we discovered 
that a composite e — m object has a meaningful characterization within the model as 
well, at least in that it has fermionic statistics. This begs the question, how many 
particle types exist in the toric code model and how can we identify them? 

We take an algebraic definition of a particle type: each type corresponds to a super- 
selection sector, which is a representation of the local operator algebra. In particular, 
we say that two particles (or composite objects) are of the same type 

a^b (4.10) 

if a can be transformed to b by some operator acting in a finite region. For example, 
in the toric code, two e-particles are equivalent to having no particles at all. 
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e 




e 







w, 



(4.11) 

by acting with an appropriate, geometrically bounded electric path operator Wi'^K 
We introduce the notation 

exe = l (4.12) 

to represent the fusion rule that two e-particles are equivalent to the vacuum sector 
1. In the toric code, there are 4 superselection sectors: 



1, e, m, and e — e x m 



(4.13) 
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with the fusion rules: 

exe = l e X m = e 

mxm~lexe = m (4-14) 
e X e ~ 1 m X e = e 

4.2 Mutual statistics implies degeneracy on the torus 



This is an argument due to [Einarsson (1990 ) . Suppose that there are at least two 



particle types, e and m with —1 mutual statistics. Let us define an operator Z acting 
on the ground state in an abstract fashion (not refering to the actual model) which 
creates an e pair, wraps one particle around the torus and annihilates the pair. In the 
toric code, this will be the path operator W^/'^'' = Yijei '^j for a loop I winding one of 
the nontrivial cycles on the torus, but we need not know that specifically. 

We can define another operator X that creates a pair of the other type m and 
winds around the other nontrivial cycle on the torus. But now a bit of geometric 
introspection reveals that the combination, 




(4.15) 

Thus, there are two non-commuting operators acting on the ground state space £, 
and we conclude dim £ > 1. In fact, there are four such operators, each of the two 
particle types can be moved around each of the two nontrivial cycles. Working out the 
commutation relations of these operators implies that dim £ = 4. 

4.3 The toric code in a field: perturbation analysis 

We now apply a magnetic field to the toric code that will realistically allow the quasi- 
particles to hop and, unfortunately, destroy its exact solvability (see |Tupitsyn et al. (2008 )) 
To wit: 

s P j 

For example, with ~ but hz ^ 0, we can view the perturbation as an electric 
path operator of length 1 on each edge. Hence, it can cause charge pair creation and 
annihilation (at an energy cost ~ 4Je) or hop existing charges by one lattice displace- 
ment, at no cost. For small hz this provides a nontrivial tight-binding dispersion to 
the charges, 

e{q) ~ 2Je — 2hz{cosqx + cosQy) (4-17) 

but does not close the gap or lead to a change in the topological degeneracy of the 
ground state in the thermodynamic limit. 

At large hz ^ Je, Jm, the model should simply align with the applied field as a 
paramagnet. Clearly, in this limit the topological degeneracy has been destroyed and 
wc have a unique spin-polarized ground state. The phase transition can be understood 
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Fig. 4.1 Numerically determined phase diagram of the toric code in a field from 



Tupitsyn et al. (20081. (I) labels the topological phase, (II) and (III) the vortex and charge 



condensates {i.e. the paramagnetic phase). The numerics were done using discrete imaginary 
time with a rather large quantization step. 

from the topological side as a bose condensation of the charges, which proliferate as 
hz increases. 

The same argument is applicable if 3> Je,Jm- If increases while hz = 0, 
then vortices condense. However, the high-field phase is just a paramagnet, so one 
can continuously rotate the field between the x- and z-direction without inducing a 
phase transition. Thus, the charge and vortex condensates are actually the same phase! 



This property was first discovered by Fradkin and Shenker (1979p for a 3D classical 
Z2 gauge Higgs model, where it appears rather mysterious. 

4.4 Robustness of the topological degeneracy 

The splitting of the ground state levels due to virtual quasiparticle tunneling is given 

by 



SE ~ Ae-^/« 
This follows from the effective Hamiltonian 

^^cff = —{tizZi + tizZi + t\x^\ + iix^-i) 
where the Z^, Xi operators are the winding loop operators of Sec 



(4.18) 



(4.19) 



Physically, this 

is simply a statement of the fact that the only way to act upon the ground state is 
to wind quasiparticles around the torus. This is a process exponentially suppressed in 
system size. 

4.5 Emergent symmetry: gauge formulation 

There are two ways to introduce symmetry operators in the perturbed toric code 
model. 
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1. One can define loop operators (e.g. Zi, Z2, Xi, X2), the definition of which depends 
on the actual ground state of the perturbed Hamiltonian. This is similar to the 
definition of the operator Y in the ID case of Sec. 12.51 which also requires detailed 
knowledge of the ground state. 

2. One can exploit gauge invariance by rewriting the model in a gauge invariant 
form. This can be done for any spin model by introducing redundancy. In this 
case, the symmetry does not depend on the model but is only manifest in the 
topological phase. 

We will take the second approach in order to avoid the difficulty of defining the 
appropriate loop operators and also to introduce the important gauge formulation of 
the model. To gauge the model we proceed in steps: 

1. Introduce one extra spin fiy per vertex that always remains in the state 

;^(lt) + U))- (4.20) 

This state is characterized by the condition 

= (4.21) 

where fx^ is the Pauli spin matrix for the spin fi at vertex v. 

2. Change spin operators from i7„„,/it, to a-uvji^v where we represent the classical 
value of each old spin Suv as m„s„t,?fit,. Here s„„ is the spin on the edge connecting 
u and V and m„, m^, are the classical values of the new spins {i.e. the labels in 
the basis). 

Thus the complete transformation is given by 

~x 

K 

K n (4.22) 

jGstar(u) 

and the constraint Eq. (|4.2ip becomes the standard Z2 gauge constramt: 

filAu\^) = \^). (4.23) 
On states satisfying the gauge constraint Eq. ()4.23p . Au = Au = /i^. Therefore, 

A.^\^>)=K\^) (4.24) 
and we can rewrite the Hamiltonian as 

H ^ - JeJ2i^v - Jrn J2 ^ J2 {hx^l^ + Kfl^KvM (4.25) 

V P (u,v) 

subject to the gauge constraint. 

Viewed as a standard Z2 gauge theory, the protected topological degeneracy of the 
ground state is physically familiar as the protected degeneracy associated with the 
choice of flux threading the 2g holes of the genus g surface. 



5 

The honeycomb lattice model 



We now investigate the properties of another exactly solvable spin model in two dimen- 
sions, the honeycomb lattice model. This model exhibits a number of gapped phases 
that are perturbatively related to the toric code of the previous sections. Moreover, in 
the presence of time-reversal symmetry breaking terms, a new topological phase arises 
with different topological properties, including nontrivial spectral Chern number. An 
extended treatment of the properties of this model with much greater detail can be 



found in Kitaev (2006) 



In the honeycomb lattice model, the degrees of freedom are spins living on the 
vertices of a honeycomb lattice with nearest neighbor interactions. The unusual feature 
of this model is that the interactions are link orientation dependent (see Fig. 15. ip . The 
Hamiltonian is 

X links y links z links 

We might expect this model to be integrable because [H, Wp] = for an extensive 
collection of plaquette operators 



w„ 



(5.2) 




Fig. 5.1 The honeycomb model has spins living on the vertices of a honeycomb lattice 
with nearest neighbor interactions that are link-orientation dependent. a;-links have a^a^ 
interactions, y-links have a^a^ interactions and z-links have a^a^ interactions. 
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where the spins and labels follow from the figure for each plaquette. Unfortunately, 
this is not quite enough: there are two spins but only one constraint per hexagon so 
that half of each spin remains unconstrained. In fact, the remaining degrees of freedom 
are Majorana operators! 

5.1 A (redundant) representation of a spin by 4 Majorana 
operators 

We consider a collection of four Majorana operators c, 6^,5^ and that act on the 
4-dimensional Fock space ^. We define the following three operators 

ay = ibyc 

ay = ibyc. (5.3) 

These operators do not obey the spin algebra relations on the full Fock space, but we 
clearly have two extra dimensions of wiggle room. In fact, the physical state space is 
identified with a two-dimensional subspace £ C 5^ given by the constraint 

Z?!*) = I*), where D = b'^byb^c (5.4) 

Within £, the tr" act as ct" act on the actual spin. Of course, ct" also act on £-^, but 
we can ignore these states by enforcing the constraint. 

To be careful, we need to check two consistency conditions: 

1. ct" preserves the subspace £, which follows from [tr", D] = 0. 

2. The a" satisfy the correct algebraic relations when restricted to £. For example, 

a'^aya' ^ {ib''c){ib»c){ib'-c) = i^{-i)b''b»b'c^ ^io^i (5.5) 

where the last equality only holds in the physical subspace £. 

5.2 Solving the Honeycomb Model using Majoranas 

We now use the Majorana representation of spins just introduced to rewrite each 
spin of the entire honeycomb model as in Fig. 15.21 This will greatly expand the 2-^^- 
dimensional Hilbcrt space to the Fock space ^ of dimension 2^^, but the physical 
space £ C 3^ is fixed by the gauge condition 

Dj\-9)^\^) for all j (5.6) 

where Dj — b^bybjC. We define a projector onto £ by 

j ^ 

In the Majorana representation, the Hamiltonian (|5.ip becomes 
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Fig. 5.2 Majorana representation of honeycomb model. 

Ajk = '^Ja(j,k)Ujk 

u,k = ibf^-'^^bf'"^ (5.8) 

where a(j, k) = x,y, z is the direction of the Unk between j and k. 

We have suggestively written the Hamiltonian H as if it were a simple quadratic 
fermion Hamiltonian as in Sec. 12.41 but of course Ajk is secretly an operator rather 
than a real skew-symmetric matrix. However, each operator enters only one term 
of the Hamiltonian and therefore Ujk commute with each other and with HI Thus, we 
can fix ujk = ±1, defining an orthogonal decomposition of the full Fock space: 

^ = 0^n, where 1^-) e;?„ iflFu.fel*) ==ii,fc|*) V j, fc (5.9) 

u 

Within each subspace we need to solve the quadratic Hamiltonian 

Hu = - ^ AjkCjCk 

Ajk = 2Ja(j,k)Ujk (5.10) 

which we know how to do in principle. On the other hand, the integrals of motion Wp 
(the hexagon operators) define a decomposition of the physical subspace £ labeled by 
the eigenvalues Wp = ±1: 

£ = 0£», where \-^) e £,.u, iS Wpl'^) ^ wp\^) y p (5.11) 

w 

We can relate these two decompositions by expressing Wp in the Majorana represen- 
tation and noting that within the physical subspace 



U,k)edp 



(5.12) 
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Thus, we find 

= Usidu (5.13) 

where Wp = Uu,k)edp '^jk- 

So we have a procedure for finding the ground state of the honeycomb modch 

1. Fix Wp = ±1 for aU p. 

2. Find Ujk satisfying 

Wp= Ujk- (5-14) 

There is a smah subtletly here in that Ujk = —u^j so we must be careful about 
ordering. We can consistently take j in the even sublattice of the honeycomb and 
k in the odd sublattice in equation (|5.14p . 

3. Solve for the ground state of the quadratic Hamiltonian ()5.10p . finding the energy 
E{w). 

4. Project the found state onto the physical subspace {i.e. symmetrize over gauge 
transformations) . 

5. Repeat for all w; pick the w that minimizes the energy. 

If there were no further structure to E{w), this would be an intractable search 
problem in the space of Wp. Fortunately, due to a theorem by Lieb (|1994p . the ground 
state has no vortices. That is, 

E{w) = min if Wp = 1 V p (5.15) 

Using this choice of Wp, it is easy to solve the model and produce the phase diagram 
of Fig. 15.31 The gapless phase has two Dirac points in the fermionic spectrum. 

5.3 Fermionic spectrum in the honeycomb lattice model 

We just need to diagonalize the Hamiltonian 

-ffu = - ^ AjkCjCk 

f +1 if 7 S even sublattice ,r ^r.\ 

"^'= = 1-1 otherwise ^^'^^^ 

This is equivalent to finding the eigenvalues and eigenvectors of the matrix iA. Since 
the honeycomb lattice has two sites per unit cell, by applying the Fourier transform 
we get a 2 X 2 matrix A{q): 



zA{q) 



ifiq)\ 
-^f{q) ) 



e{q)^±\f{q)\ (5.17) 

where f{q) is some complex function that depends on the couplings J^, Jy, Jz. In the 
gapless phase (phase B in Fig. l5.3p . f{q) has two zeros which correspond to Dirac points 
(see Fig. 15. 4[) . At the transition to phase A, the Dirac points merge and disappear. 
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J= (0,0,1) 




/= (1,0,0) /= (0,1,0) 

Fig. 5.3 Phase diagram of honeycomb model. This is a slice through the positive octant in 
J coupling space along the Jx + Jy + Jz ~ 1 plane. The other octants are analogous. 




Fig. 5.4 Direct and reciprocal lattices of the honeycomb. The points ±(f* are the two Dirac 
points of the gapless phase B. 

5.4 Quasiparticle statistics in the gapped phase 

It appears that there are two particle types: fermions and vortices (hexagons with 
Wp = —1). The vortices are associated with a Z2 gauge field, where Uj^ plays the role 
of vector potential. Taking a fermion around a vortex results in the multiplication of 
the state by —1 (compared to the no- vortex case). However, the details such as the 
fusion rules are not obvious. 

Let us look at the model from a different perspective. If Jx = Jy = 0, J, > 0, 
the system is just a set of dimers (see Fig. 15. 5|) . Each dimer can be in two states: ft 
aird ii- The other two states have 2 higher energy. Thus, the ground state is highly 
degenerate. 

If Jx, Jy ^ Jz, we can use perturbation theory relative to the noninteracting dimer 
point. Let us characterize each dimer by an effective spin: 

l^) = ITT); l^) = IU). (5.18) 

At 4th order of perturbation theory, we get: 

const -^^g, (5.19) 

2 P 
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Fig. 5.5 The vertical dimers on tlie iioneycomb lattice themselves form the edges of a 
(dashed) square lattice. The plaquettes of alternate rows of the hexagonal lattice correspond 
to the stars and plaquettes of the square lattice. This is weak breaking of translational sym- 
metry. 

where p runs over the square plaquettes of the dimer lattice (see Fig. 15. 5[) and 

= (5-20) 

is a plaquette operator on the effective spin space | tl" ), | 4 )• By adjusting the unit 
cell and rotating the spins, wc can reduce this Hamiltonian to the toric code! 

The vertices and plaquettes of the new lattice correspond to alternating rows of 
hexagons. Thus, vortices on even rows belong to one superselection sector and vortices 
on odd rows to the other. It is impossible to move a vortex from an even row to an 
odd row by a local operator without producing other particles {e.g. fermions). The 
fermions and e — m pairs belong to the same superselection sector, e, though these are 
different physical states. 



5.5 Nonabelian phase 

In the gapless phase B, vortex statistics are not well-defined. However, a gap can 
be opened by applying a perturbation that breaks the time-reversal symmetry, such 
as a magnetic field. Unfortunately the honeycomb model in a field is not exactly 



solvable. Yao and Kivelson (2007) studied an exactly solvable spin model where the 
time-reversal symmetry is spontaneously broken, but wc will satisfy ourselves by in- 
troducing a T-breaking next nearest neighbor interaction on the fermionic level (which 
can be represented by a 3-spin interaction in the original spin language). 
Written in terms of Majorana fermions, we consider the Hamiltonian 

^ = i E ^J'^^^C'^- (5.21) 



where Ajk now has chiral terms connecting Majoranas beyond nearest neighbor in the 
honeycomb lattice (see Fig. I5.6p . After Fourier transforming, we find 
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Fig. 5.6 Picture of chiral interaction matrix Aj^- Forward arrows correspond to positive 
entries in the skew-symmetric real matrix Ajk- Solid arrows are the interactions of the original 
honeycomb model; dashed arrows give the time-reversal symmetry breaking perturbation. 



with the massive dispersion relation 



(5.22) 



(5.23) 



Within this massive phase, we wiU find nontrivial topological invariants of the 
quasiparticlc spectrum. Let iA{q) be a nondcgencratc Hermitian matrix that conti- 
nously depends on q. In our case, A acts in C^, but in general it can be C" for any 
n. Let us keep track of the "negative eigenspace" of iA{q): the subspace £{q) C C" 
spanned by eigenvectors corresponding to negative eigenvalues. For matrix (|5.22p . 
dim £{q) = 1. This defines a map F from momentum space (the torus) to the set of 
m-dimensional subspaces in C". More formally: 



F : r 



U{n)/U{m) X U{n-m) 



(5.24) 



This map F may have nontrivial topology. 

In the honeycomb model with T-breaking, we have n = 2, m = 1 and U{2)/U{\) x 
f7(l) = = 5*2 is the unit sphere. Thus, F : — > S"^ and for the matrix iA{q) 
of Eq. (j5.22|l . F has degree 1. That is, the torus wraps around the sphere once. More 
abstractly, 2,{q) defines a complex vector bundle over the momentum space T^. This 
has an invariant Chern number which in this case is v ^ \. 

What is the significance of the spectral Chern number? It is known to characterize 
the integer quantum Hall effect, where it is known as the "TKNN invariant". For a 
Majorana system, there is no Hall effect since particles are not conserved. Rather, the 
spectral Chern number determines the number of chiral modes at the edge: 



V = (^ of left-movers) — (# of right-movers). 



(5.25) 



5.6 Robustness of chiral modes 



A chiral edge mode may be described by its Hamiltonian: 
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Hedge = ^ J fl{x)dxfi{x)dx 

where f}{x) is a real fermionic field. That is, 

f}{x)f]{y) + fi{y)fi{x) = 26{x - y). 
At temperature T, each mode carries energy current 

24 



(5.26) 



(5.27) 



(5.28) 



The easiest explanation of this is a straightforward 1-D fcrmi gas calculation: 



ede 



2tt 



24 



1 + e^/T 



(5.29) 



However, it is useful to reexamine this current using conformal field theory (CFT), 
in order to understand better why the chiral modes are robust. We consider a disc 
of B phase extended into imaginary time at temperature T. That is, we have a solid 
cylinder with top and bottom identified: 




We have obtained a solid torus whose surface is a usual torus. The partition function 
is mostly determined by the surface. 

Let the spatial dimensions be much greater than ^ . From this point of view, the 
cylinder looks more like: 




According to the usual CFT arguments, we have 
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T = i h twist (5.30) 

V 

Twisting the torus changes the partition function by: 

T i—i- T + 1 

Z^Ze^""'^. (5.31) 

On the other hand, the twist parameter (Rer) couples to the some component of the 
energy-momentum tensor, namely, T^.^, which corresponds to the energy flow. This 
relation implies that 

I=^{c~c)T\ (5.32) 

The chiral central charge, c — c, does not depend on the boundary conditions. Indeed, 
the energy current on the edge cannot change because the energy cannot go into the 
bulk. 
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